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Abstract 

The -Bs —7- decay plays an outstanding role in tests of the Standard Model 

and physics beyond it. The LHCb collaboration has recently reported the first 
evidence for this decay at the 3.5 a level, with a branching ratio in the ballpark of 
the Standard Model prediction. Thanks to the recently established sizable decay 
width difference of the Bs system, another observable, A^^^, is available, which 
can be extracted from the time-dependent untagged Bg — ^ rate. If tagging 

information is available, a CP-violating asymmetry, S^^, can also be determined. 
These two observables exhibit sensitivity to New Physics that is complementary 
to the branching ratio. We define and analyse scenarios in which these quantities 
allow us to discriminate between model-independent effective operators and their 
CP-violating phases. In this context we classify a selection of popular New Physics 
models into the considered scenarios. Furthermore, we consider specific models 
with tree-level FCNCs mediated by a heavy neutral gauge boson, pseudoscalar 
or scalar, finding striking differences in the predictions of these scenarios for the 
observables considered and the correlations among them. We update the Standard 
Model prediction for the time-integrated branching ratio taking the subtle decay 
width difference effects into account. We find (3.56 ±0.18) x 10~^, and discuss the 
error budget. 
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1 Introduction 



The rare decay Bg — t- is very strongly suppressed within the Standard Model (SM). 

As this need not be the case for SM extensions, which could dramatically enhance it, 
it has already for decades played an important role in constraining such extensions and 
giving hope for seeing a clear signal of New Physics (NP) Using the parametric 
formula for its branching ratio from Ref. J2], and updating the values for the iJ^-decay 
constant Fb^ [3j and the _Bs-lifetime tb^ [4]7 we find 

BR{Bs /i+/i")sM = (3.25 ± 0.17) x 10"^ (1) 

Thus, within the SM, only about one in every 300 million mesons is predicted to 
decay to a pair of muons. 

Concerning the measurement of this branching ratio, a complication arises due to the 
presence of B^-B^ oscillations [s] . In particular, LHCb has recently established a sizable 
value of the decay width difference AF^ between the Bg mass eigenstates (gI: 



ys = ^, ^ = = 0.087 ± 0.014, (2) 



L ^ ^ H 

-1 



where = r^^ denotes the average Bg decay width. This quantity enters the time- 
integrated decay rate, which is at the origin of the measurement of the "experimental" 
branching ratio BR(i?s — /i^/i^). It is related to the "theoretical" branching ratio, 
referring to ?/s = 0, as follows Isl: 



BR{Bg 



1 + ^^yVs _ 



BR(5, ->/iV~), (3) 



where the observable A'^^' eQuals +1 in the SM. Using the numerical value for the SM 
branching ratio in ([T]) and the experimental value for yg in ^ gives 

BR{Bg ;U+/i")sM = (3.56 ± 0.18) x 10"^ (4) 

which is the reference value for comparing the time-integrated experimental branching 
ratio with the SM. 

Over the last decade we have seen the upper bounds for the Bg — )■ ^+fi^ branching 
ratio continuously move down thanks to the CDF and DO collaborations at the Tevatron 
and the ATLAS, CMS and LHCb experiments at the LHC (for a review, see Ref. [7]). In 
November 2012, the LHCb collaboration reported the first evidence for the Bg — )■ fi'^fi~ 
decay at the 3.5o" level, with the following branching ratio Isl: 



BR{Bg /iV")LHCb = (3.2+i:^) X 10^' G [1.1, 6.4] x 10"" (95% C.L.). (5) 

The agreement with Q is remarkable, although the rather large experimental error still 
allows for sizable NP contributions. It is already obvious at this stage, however, that 
this will be a challenging endeavour. 



^This should be compared with BR{Bs fJ-+^J.~)sM = (3.23 ±0.27) x 10'^ in Ref. [|, implying that 
the central value remains practically unchanged but the error has decreased significantly. We discuss 



the reason for this change in Subsection 2.3 
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As emphasized in Ref. [5], and discussed and illustrated in detail below, the observable 
entering ([s]) is sensitive to NP contributions. It is thus a complementary observable 
to the Bg —7- ^^fi~ branching ratio, offering independent information on the short- 
distance structure of this decay. The observable can be extracted from the untagged 
data sample, for which no distinction is made between initially present or B^ mesons, 
once enough decay events are available for the decay-time information to be accurately 
taken into account |5|[9]. The conversion factor in ([s]) would then be determined from 
the data, thereby also allowing the extraction of the theoretical Bg — )■ fi~^fi~ branching 
ratio. 

If tagging information is included, requiring even more events to compensate the 
efficiency of distinguishing between initially present B^ or B^ mesons, a CP-violating, 
time-dependent rate asymmetry can be measured. As we assume that the muon helicity 
will not be measured, this rate asymmetry is governed by another, third observable, 
iS^^, which vanishes in the SM but is very sensitive to new CP- violating effects entering 
Bg — 7- /i''"/^" in extensions of the SM [s]. Analyses of CP- violating effects in Bs(d) 

decays, which neglected AF^ effects, were performed for various models of NP in 



Refs. 10 13 . An analysis of Z' models that including AF^ effects was performed in 
Ref. ^ 

The observables A'^^ ^f^f^ both depend on the CP-violating B'^-Bg mixing phase 

^ _ ^SM ^ ^NP ^ _2A2r/ + 0f , (6) 



s 



where the numerical value of the SM piece, involving the Wolfenstein parameters A and 
r] of the CKM matrix, is given by -(2.08 ± 0.09)°. The LHCb analysis of CP-violation 
in the Bg — )■ J/ip(f) decay currently gives the most precise experimental determination of 
this phase [g]: |^ 

0^ = - (0.06 ± 5.99)° . (7) 

Hadronic uncertainties from doubly Cabibbo-suppressed penguin effects have been ne- 
glected in this measurement; these corrections have to be controlled once the experimen- 



tal precision improves further 16 



As the measurements of the observables A'^ and refer to the era of the LHCb 
upgrade, we can assume that (j)g will be known precisely once data for these observables 
become available. The goal of the present paper is to investigate and illustrate how the 
experimental knowledge of the trio from Bg — t- fj^^fi~, 

BR(5, ^/iV), ^Ar, (8) 

will shed light on the possible presence of NP in this decay, which cannot be obtained 
on the basis of the information on BR(i?s — t- fj^~^fi~) alone. 

Our paper is organized as follows: in Section [2] we recall the definitions of the ob- 
servables in ([s]) and discuss their various properties. In Section [s] we introduce various 
scenarios for NP, classifying them in terms of four general parameters which can in prin- 
ciple be calculated in any fundamental model and are directly related to the physics of 



■^We avoid averages for 0^ that include decays such as — J /ipf o{980) because of the unsettled 



hadronic structure of the /o(980) state 15 
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Bs — )■ jji'^ jj," . In this context we classify a selection of popular NP models into the con- 
sidered scenarios. In Section |4] we consider three classes of specific NP models. The first 
one with tree level neutral gauge boson contributions to FCNC processes that could rep- 
resent Z' models or models with flavour violating Z couplings. The second one in which 
the role of gauge bosons is taken over by scalar or pseudoscalar tree-level exchanges. Fi- 
nally, we present a third model in which both a scalar and a pseudoscalar with the same 
mass couple equally to quarks and leptons. We demonstrate how the measurements of 
the observables in Q can distinguish between these three classes of models. In Section |5] 
we summarize the highlights of our paper. 



2 Observables of Bg fJ^^fJ' 



2.1 Basic Effective Hamiltonian 

The general model-independent low-energy effective Hamiltonian for a. Bg £^i^ decay 



IS 



17 19 



w,s,p •\ 



HefT = <; vzvt, ^ {a + c[ o[) + h.c ;> , (9) 

where the operators are 

Os = mt{sPRb){ll), O's = mh{sPLb){ll), 

Op = m,{sPnb){h5l), O'p = m.isPMh^l) (10) 

and a = e'^/An is the QED fine structure constant. The observables that we will calculate 
below can each be expressed in terms of the following combinations of Wilson Coefficients: 



safe) 1^1^'- 

In the SM Cio = Cf^ and C[q, and &!) are all negligibly small, so that P^^ = 1 
and = 0. The Wilson coefficient Cio is given in the SM as follows 

CiT = -r]Y sin-2 OwYoixt) = -4.134, (12) 

where Yo{xt) is a one-loop function with xt = mf/M^ [20), and the coefficient rjY is a 



QCD factor that for rrit = mt{mt) is close to unity: rjy = 1.012 21,22 . The scalar 
Wilson coefficients are thereby related to the parameter S by the numerical factor 



mk{Cs - C's) = -0.130 X S. (13) 

lile Cio and C(q a 

sion GeV~^. 



Note that while Cio and C(q are dimensionless, the coefficients C^'^ and Cp"* have dimen 
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2.2 Time-Dependent Rates 

The time- dependent rate for a 5° meson decaying to two muons with a specific hehcity 
\ = L, R is given by 



Cl cos(AM, t) + S^^ cos(AM, t) 



m 



\P\' + \S\') 



Bs 



+ cosh( ^) +^^'^rsinh( ^ 



X e 



(14) 



where r^^ = 2/{Tii + Fl) is the Bg mean hfetime and t/^ is defined in (j2]). 

The time- dependent rate for a 5° meson is obtained from the above expression by 
replacing C^^ — )■ — and — )■ — The time-dependent observables for both rates 



can be expressed in terms of the parameters defined in (11) as [5 ,23 
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(15) 
(16) 
(17) 



The phase represents the CP-violating NP contributions to B^-B^ mixing. It in 



fluences the mixing-induced CP asymmetry in the B^{Bg) — )■ ^/'^ decays 
latter given by 
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• sin 0, = sin(2 1/3,1 



Vts = -\Vts\e~ 



with the 



(18) 



with (3s — —1° . In Ref. hi and the papers reviewed there 0, ^ = 2^9^^ 

Only the observable C^^ is dependent on the helicity of the final state i.e. it depends 

on the parameter rjx = {+1: L; —1: R}. The presence of the observable A'^/' is a 
consequence of the sizable Bg decay width difference AP,,. 

In practice the muon helicities A are very challenging to measure. If no attempt is 
made to disentangle them, then we measure their sum: 



p(s°(t) ^ /iV") = Yl ^^^sit) ^ fitf^x), 

X=L,R 

p(s°(t) /iV-) = E r(^°(t) ^ ^^tf^-,). 



(19) 



X=L,R 



helicity, cancels in both sums l5 



Observe from equations (|14|) and (|T5|) that C^^, which was dependent on the muon 



The Bg — !■ n~^fJ^ helicity-summed time-dependent untagged rate is then given by 

M,l sin^ 



'-'10 ^tsytbl i'Ss^Bs^ 



Ami 



47r5 

X (|pp+i^n 



X e 



[cosh t/ Tb, ) + ^^'^^ sinh {y^ t/ J] . 



Similarly, the helicity-summed time-dependent tagged rate asymmetry is 
r(50(t) ^ - r(50(t) ^ _ sin(AM,t) 



r(50(t) ^ + r(50(t) ^ cosh(?/,t/rBj + AZ sinh(y,t/rB, 



(20) 



(21) 



It is important to clarify that although there is no explicit term for direct CP violation 
in the rate asymmetry, this does not mean that the absolute values squared of S^^ and 
A^^Y necessarily sum to one. These two observables also have an implicit dependence on 
C^^, the rate asymmetry for 5° and i?^ decays to the specific helicity muon final states. 
This gives the relation 



MM 
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1 - \c 



X |2 
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(22) 



Thus if there are no new CP-violating phases in the mixing or decay amplitudes, ifp = 
ips = = such that = 0, ATr does not have to take its SM value of 1. The 



presence of a non-negligible scalar operator Og\ so that l^l ^ 0, is sufficient to ensure 



that A^Y 7^ "^^^ ^1^° be seen from (17) 



The dependence on Fb^ cancels in both A'^^ ^^"^ '^mm' consequently they are 
theoretically clean. Moreover, these observables are also not affected by the ratio fd/fs 
of fragmentation functions, which are the major limitation of the precision of the Bg — > 
branching ratio measurement at hadron colliders 25 . As A'^y does not rely on 
flavour tagging, which is difficult for a rare decay, it will be easier to determine. Given 
enough statistics, a full flt to the time-dependent untagged rate will give A'^^- With 
limited statistics, an effective lifetime measurement may be easier, which corresponds to 
fltting a single exponential to this rate. For a maximal likelihood fit, the Bg IJ''^fJ'~ 



effective lifetime is equal to the time expectation value of (20) [9]: 

_/,°°t(r(i?,(t)^^V-))rft 



T, 



/,°°(r(i?,(t)^/i+/i-))rft 

The untagged observable is then given by 

il-y')T„„-a 



(23) 



/IMM 



AF 
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(24) 
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2.3 The Branching Ratio 



A Bg fi'^fi branching ratio measurement amounts to counting all events over all 
(accessible) time, and is thus defined as the time integral of the untagged rate given in 

m mm 



{T{B,{t) ^ fi+fi-)) dt. 



(25) 



LHCb has recently presented the first measurement of the Bg — )■ ^~ time-integrated 
rate fs] that we have given in (|5]). In contrast, the SM prediction for the Bg — )■ 
branching ratio in ([T]) is computed theoretically for one instant in time, namely at t = 
i.e. it neglects the effects of B^-B^ mixing. Specifically, it is given by 



jSM 



(r(i?,(t)^/iV)) 



t=0, P=l,5=0 



TB, G^M^sin^^H. 
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(26) 



an updated numerical estimate is given in ([T|. 

It is now straightforward to derive the expression 



|P|' + 1^1'. 



(27) 



BR{Bs )sM 

However, as not the theoretical but the experimental branching ratio is measured, it is 
useful to introduce the following ratio js]: 



R 



BR{Bs ii+fi-) 
BR{Bs fi+fi-)sM 
l + Us cos(2v9p - (f)^^] 



l+Vs 



\p\ 



xi\P\' + \S\') 
Vs cos{2(ps - (j)^^) 



(28) 



where the sizable decay width difference AF^ enters. The parameter R is related to R 
defined in Ref. js] by i? = (1 — ys)R- Combining the theoretical SM prediction in ([T]) 
with the experimental result in ([s]) gives 

R = 0.901°:^^ g [0.30, 1.80] (95% C.L). (29) 

This range should be compared with the SM value Rsm = 1- 

Finally we would like to explain the origin of the reduced error in To this end 
we return to the basic parametric formula (18) for the theoretical branching ratio in 
Ref. [2]. It turns out that the changes of the input parameters over the last six months 
have practically no impact on the central value obtained there. Indeed updating the 
central values of Fb^ and tb^ , we can cast this formula into the following expression: 



BR{Bg /i+/i")sM = 3.25 X 10" 



173.2 GeV 



3.07 



225MeVy Vl-500psy 



0.0405 
(30) 
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1 .5% 




0.8% 




4.0% 



Figure 1: Error budgets for the two branching ratio calculations of Bg — )■ in the 



Standard Model given in (30) (left) and (33) (right). 



(31) 



The most recent world averages for Fb^ [3] and tb^ (4] are 

Fb, = (225 ± 3) MeV, tb, = 1.503(10) ps 

to be compared with Fg, = (227 ± 8) MeV and tb, = 1.466(30) ps used in Ref. |2|. 
While the change in tb, is an experimental improvement, confirmation of the impressive 
accuracy on Fb^ is eagerly awaited. In Ref. |2] a more conservative approach has been 
used, but here we follow Ref. j3], updating also tb^- With unchanged input on Mt and Vts 
with respect to Ref. [2] we arrive at ([TJ and consequently, after including the correction 
from Ar„ at (Q. 

Now as stressed and analysed in |2}[26] additional modifications could come from 
complete NLO electroweak corrections, which should be available soon (M. Gorbahn, 
private communication). They are expected to be at the level of at most a few percent. 
Taking at face value the present error on Fb^ , the current error budget for the branching 
ratio is as follows: 



Mt: 1.5%, 



Fb, : 2.7%, TB, : 0.7%, 



4%, 



(32) 



It is also depicted in the left panel of Figure [T] Not knowing yet the correction from 
NLO electroweak effects, evidently the error on |V^^^ts| is now the largest uncertainty 
but this assumes that the error on Fb^ is indeed as small as obtained in Ref. [3]. 

While the small error on Fb^ is expected to be consolidated soon, the decrease of the 
error in \Vts\ appears to be much harder. In this context it should be recalled that the 
branching ratio in question can also be calculated by using the mass difference AM, 
The updated parametric formula (13) of the latter paper reads 



27 



BR{Bs ;U+/i")sM = 3.38 x 10" 



173.2 GeV 



1.500psy \B 



1.33 



AM. 



17.72/ps) ■ 
(33) 

We note that among the uncertainties in (32), the largest two are absent and the uncer- 
tainty due to Mt is reduced to 0.8%. The uncertainty due to AM^ is negligible [i]. The 
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error budget for this expression reads 



Mt : 0.8%, TB^ : 0.7%, AM, : 0.3%, Bb, : 4.5%, (34) 
where we used Bbs = 1.33(6) [28j. It is pictured in the right panel of Figure [ij The 



latter uncertainty is expected to be reduced significantly in the coming years so that (33) 
could remain to be the most accurate estimate of the branching ratio in question under 
the assumption that there are no NP contributions to AM, and the SM reproduces its 
experimental value. In view of this tacit assumption in (33), we prefer to use ([T| as the 
present best estimate of the theoretical branching ratio. On the other hand, by using 



(33) we find, after the inclusion of AF, the correction, 

ER{Bs /U+/i~)sM = (3.71 ± 0.17) x 10"^ (AM,) (35) 

which agrees very well with (0). This updates the estimate (3.2 ±0.2) x 10"^ in Ref. 
where AF, effects where not included and other input, in particular the value of tb^, was 
different than now. 



3 Constrained Scenarios and Their Phenomenology 
3.1 Preliminaries 

Experiments have started honing in on the B^ — )■ fi^fj^~ time-integrated rate, or branching 
ratio, for which the observable R parameterises possible NP contributions. Next in 
line is a time-dependent analysis, first without tagging, giving ^^p, and then with 
tagging, giving iS^^. The end result will be three experimental observables, which, if 
there are scalar operators contributing to the decay mode, can each contain independent 



information (see the discussion around (22)). 

With the phase (p^^ already significantly constrained by the current data ([t]), these 
three observables depend on four unknowns: 

|P|, (fp, \S\, LPs- (36) 

Therefore we cannot in general solve for all of these model-independent NP parameters 
by considering the decay Bg — )■ fi'^fJ'~ alone. One solution is to invoke other b — )■ 
sfi~^fi~ transitions like the decays B — t- Kfi'^fi~ and B — )■ K*fi^fi~. In particular, as 



emphasized in Ref. 30 , observables in i? — )• Kfi'^fi~ are sensitive to Cs,p + C'g p, rather 
then differences of tEese coefficients, thereby allowing additional complementary tests 
and in principle the determination of all Wilson coefficients. But present form factor 
uncertainties in these decays do not yet provide significant new constraints on scalar 
operators relatively to the ones obtained from Bg — )■ fj,'^fj.~. In any case such analysis 
would be beyond the scope of the present paper. 



In the spirit of the analysis of Z' contributions to FCNC processes in Ref. 14 , where 
various scenarios for Z' couplings to quarks have been considered, and an analogous 
analysis for tree-level scalar and pseudoscalar contributions to FCNC processes |31j, we 
will consider various scenarios for S and P that will allow us to reduce the number of free 
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NP parameters and eventually, with the help of future data, uniquely determine them. 
Our scenarios are motivated by generic features of NP models and, as we will see below, 
they result in a distinct phenomenology for the observables i?, v4^p and 5^^. In the 
present section our analysis is dominantly phenomenological, although we discuss the 
motivation behind each scenario and the characteristic features of its phenomenology. 
Moreover we indicate what kind of fundamental physics could be at the basis of each 
scenario considered and we survey specific models of NP and categorise them into the 
scenarios that we will list now. 

The five scenarios to be considered are as follows: 

(A) S = 

(B) P = 1 

(C) P±S = 1 

(D) ipp,ifis e {0,7r} 

(E) P = 0. 

The scenarios are intended to be limiting cases, i.e. although we are not aware of a 
model that predicts P = 0, P ~ is conceivable and the resulting phenomenology will 
be approximately the same. 



3.2 Scenario A: S = 



3.2.1 General Formulation 

This scenario is realised if Cs — C'„ 



s 



0, 



leaving cfj 



and Cp"* free to take non-SM 



values as well as CP-violating phases. Thus models with only new gauge bosons or 
pseudoscalars naturally fall into this category and consequently, as we will see below, 
this scenario includes a number of popular BSM models. Also models with scalars can 
quahfy, provided the scalars couple left-right symmetrically to quarks so that Cs = C'g. 

In this scenario the rate asymmetry between and P° decays to the individual 
muon helicities vanishes: C^^ = 0. Therefore the two time-dependent observables do not 
carry independent information, being bound by the constraint 



'MM I 



+ 1^: 



Arl 



1. 



sin(2(^/ 



Specifically, 

while the branching ratio observable is given by 

"l + l/,cos(2y;p-0NP) 



R= P 



(37) 



(38) 



(39) 



Formulae (38) and (39) are the basic expressions for this scenario. The three observ- 



ables in (38) and (39) are given in terms of two unknowns: |P| and (pp. As we assume 
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knowledge of and S^^ will allow an unambiguous extraction of the phase 2ipp, 

which in turn, with the help of R, will provide the value of \P\. 

The P parameter can also conveniently be expressed as P = 1 + P with 

P = iPle- = - + < ( \ (40) 

where 

5C,o ^ C,o - Cf^. (41) 

In this notation all NP effects are contained in the parameter P. In the left panel of 
Figure [2] we show the correlations between R and Scenario A using this notation. 

We have varied P G [0,1], and most importantly show the strong dependence on the 
phase (fip. As will be discussed in detail in Section |4| the requirement for new gauge 
bosons or pseudoscalars to satisfy the Bg mixing constraints implies that (fp ~ it/ 2 or 
ifp ~ 0,7r, respectively. 

Note that in the case of no new phases, ipp G {0, tt}, and = 0, 

^Ar = l, 5^M = 0, P=|P|^ (42) 
While the first two results coincide with the SM, NP effects can still arise in R. 

3.2.2 Examples of Models 

Constrained Minimal Flavour Violation (CMFV) 

In the CMFV scenario it is assumed that new low-energy effective operators beyond 
those present in the SM are very strongly suppressed and that flavour violation and 
CP- violation are governed by the CKM matrix |32^|33j. Thus all the Wilson coefficients 
aside from Cio are zero, and Cio is real. This translates into Scenario A, with the added 



restrictions that ifp = = 0. Consequently the formula (42) applies and NP enters 
only through the ratio R. 

Littlest Higgs Model with T-Parity (LHT) 

Similar to CMFV, only SM operators are relevant in this framework but due to the 
presence of new phases in the interactions of SM quarks with mirror quarks, CP asym- 



metries can differ from the SM ones. Therefore general formulae (38) and (39) apply 
here. Typically BR(Ps is predicted to be larger than its SM value but it can 

only be enhanced by 30% at most |34|. A significant part of this enhancement comes 
from the T-even sector that corresponds to the CMFV part of this model, while \S^(j)\ 
and \Sfj_f^\, governed by new phases in the mirror quark sector, are at most 0.2. 

Z' Models and RSc 



As demonstrated in Ref. 14 , larger effects than in LHT can be found in Z' models 
with tree-level FCNC couplings. If new heavy neutral gauge bosons dominate NP con- 
tributions to FCNCs, S" = in these models, placing them automatically in Scenario 



10 



A. In contrast to CMFV and LHT, the presence of new operators implies a rather rich 
phenomenology. Yet, in the absence of S the time-dependent observables are only sen- 



sitive to new CP violating phases and are not independent of one another. In Ref. 14 



correlations between the S^^ observable and AF = 2 observables have been found for 
different scenarios for Z' couplings with the size of effects that could be measured in the 
future provided the masses of these new gauge bosons do not exceed 2-3 TeV. We will 
return to this scenario in Section |4] showing results complementary to the ones presented 



in Ref. 14 



Smaller, but still measurable, effects have been found in 331 models in which new 
CP phases are present but no new operators |35|. Here NP effects in Bg — )■ fi~^fi~ are 
comparable to the ones in the LHT model provided the mass of the new neutral gauge 
boson does not exceed 2 TeV. 

Finally we mention the Randall-Sundrum model with custodial protection in which 
NP contributions to Bg — )■ fi^fi" are governed by right-handed flavour-violating Z cou- 
plings to quarks but the resulting branching ratio is SM-like, with departures from SM 
prediction at most of order 15% 136]. Larger effects are found if the custodial protection 



is absent and then left-handed couplings dominate 37 . Recently detailed analyses of 



Z couplings in similar scenarios related to partial compositeness have been presented in 



Refs. 38 40 . While having different goals than in Ref. 14 , they also demonstrate the 



power of Bg — ?■ n n in distinguishing between various NP scenarios. 



Four Generation Models 

In spite of the fact that the existence of a fourth generation seems to be very unlikely in 
view of the LHC data, in particular Higgs branching ratios, we just mention that it also 
belongs to Scenario A. NP effects in Bg — )■ fi~^fi~ can still be sizable in these models. 



See Ref. 41 and references therein. 



Pseudoscalar Dominance 

Also a model with NP dominated by tree-level FCNC contributions of a pseudoscalar 



belongs to this class. It has been analysed recently in Ref. 31 and we will present 



complementary implications of this model particularly suited to our paper in Section |4j 



3.3 Scenario B: P = 1 
3.3.1 General Formulation 

The simplest realisation of this scenario is = Cf^ and C[o = dp> = 0. However, 
pseudoscalars that couple left-right symmetrically to quarks, so that Cp = C'p, or a 
conspiracy of the form Cio — CJq = Cf^ are also allowed. The point is that in this 
scenario only scalar operators Og^ drive new physics effects in Bg — )■ In this sense 

this case is complementary to Scenario A. 

As there are scalar operators present, there is a rate asymmetry in the B^ and B^ 
decays to the individual muon helicities. Therefore the two time-dependent observables 
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P=l + \P\ e'"', s = n 
\P\<1, •fp€[l^] 
\P\<1, ¥Jpe[0,flU(3f,7rl 



|S|, 495 free; P = l 
Excluded at 95% C.L. 



1.8 



Figure 2: The correlation between the R and observables in Scenario A (left panel) 
and Scenario B (right panel). In Scenario A we have set P = 1 + P and 5 = with P 
free to vary. In Scenario B P = 1 and S is free to vary. 



do carry independent information. In this scenario the observables are given by 

COS0fP- |5|2cOs(2(/^5-0f^) 



Ar 



5, 



sin0^P_ |^|2sin(2y,^_0NP^ 



fin 



1 + 2/. 



l-y,cos(2ypg-0fP) 



(43) 



Again, with precise value of to be determined first, these three observables are 
in principle sufficient to determine the two NP unknowns, 2ips and IS*]. Consequently 
the untagged observables R and are already sufficient to determine 2{ps and IS*]. 
Moreover, if all three observables are considered, correlations between them will result 
that depend on the precise value of cf)^^ 



31 



In the right panel of Figure [2| we show the correlation between R and A'^y' ^'^^ different 
values of S* Is' 



An interesting feature is that for no CP violating phase, (ps = {0, tt}, an 
— )■ 0. But within current experimental bounds we have the 
prediction that ^^p cannot take a negative value. Moreover in this scenario IS"! < 0.5 is 
favoured. 



increase of 1 5*1 pushes A'Tr — ?■ 0. 



3.3.2 Examples of Models 
Scalar Dominance 

A model with NP dominated by tree-level FCNC contributions of a scalar belongs to 
this class. It has been analysed recently in Ref. 31 and we will present complementary 



implications of this model particularly suited to our paper in Section |4j 
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3.4 Scenario C: P ± S = 



1 



3.4.1 General Formulation 



The meaning of this scenario is clearer if we let P = 1 + P , with P defined in (40). 
Then the condition P ± S* = 1 is equivalent to P = i.e. in this scenario NP effects 
to S and P are on the same footing. If we neglect contributions to Cfg 



respect to tub^ 
Letting P = 



this scenario is realised if Cg"^ 



±C 



(') 
p ■ 



and with 



-kS for K 



5, 



±1, the time dependent observables are 

_ cos 0f P - 2 k\S\ cos((^5 - 

1-2k\S\cos(Ps + 2\S\^ ' 

-sin0NP_2fi:|^|sin(v;s-0D 



1 — 2k\S\ cos (fs + ^ISl"^ 
which are in general independent. The branching ratio observable is 
1 — 2k\S\ cosips + Sl^p + ys[cos(j)^^ — 2k\S\ cos{ips — 



R 



l + Vs 



(44) 



(45) 



In the presence of a precise value of 0^^, these three observables are sufficient to deter- 
mine the two NP unknowns ips and IS*]. Moreover, correlations between the involved 
observables characteristic for this scenario and additional tests are possible. 

The observable R is minimised by ^crit = '^(1 + 2/s)/2 and = 0, giving the lower 
bound 

R > i^. (46) 

This lower bound, without the ys and phase considerations, was first observed in Ref. [42] . 
A branching ratio measurement below this bound would thereby rule out this scenario. 

If we assume the new physics phase (p^^ in Bg mixing is known, then the purely 
untagged observables A'^^ ^ solve for S and ips- Setting = for simplicity, 
we have the expressions 



1^1 



IP 



'R{l + ys){l 



COS (/9 s = — fi;cos(<^p) 



2{1 + y,AZ) 



2P(l + l/s)(l 



P(l + l/,)^^^r 
l + VsAZ 



(47) 



iNP 



In the left panel of Figure |3j we show the correlation between R and A'^y limit 
0. Observe the lower bound on R specified in (46). If, furthermore, ipp 



{0, vr} we observe that ^4^^ can help to resolve the two possible solutions for 5* coming 
from a branching ratio measurement R. 

In the right panel of Figure [s] we show the correlation between and IS"]. Observe 
that the current measurement of R still allows a large range for both NP parameters. 
If A^^Y were measured with a negative sign it would indicate large contributions from 
NP. Moreover in this case the A^^^ sharply cuts the R contour, so that a measurement 
of A^^Y would distinguish between the magnitude and the phase of S up to the twofold 
ambiguity in 995. 
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Figure 3: Scenario C: P ± S' = 1. Left panel: the correlation between the R and „4^p 
observables. Right panel: correlation between the 15*1 = |P — 1| andips = ifp+{l + K)7i/2 
NP parameters (see text). 



3.4.2 Examples of Models 

Two Higgs-Doublet Models (2HDM), MSSM 



A 2HDM in the decoupling regime, such that M^o — Mao — Mh± ^ Mh 43 , has the 
generic feature that 

Cs = -Cp, C's = C'p. (48) 

If the couplings of the heavy Higgs bosons are not left-right symmetric, so that either 
Cs,p or C'g p are dominantrl this corresponds to Scenario C. Thus the branching ratio has 
a lower bound and a significant scalar NP contribution is indicated by negative values of 
A^^Y- ^ precise measurement of the untagged observable A^^^ distinguish the phase 
and magnitude of the NP Wilson coefficients. We will analyse a similar scenario in more 
detail in Section HI 

The above is true also for the MSSM, provided that NP contributions to vector-axial 
operators, C[q, are negligible. The MSSM has the added advantage that large tan/3 
effects, which are one way to realise the decoupling regime, can give a significant boost 



to the scalar operators 44 -46 



If the 2HDM is not in a decoupling regime, then either the physical scalar H or 
pseudoscalar may be considerably lighter than the other. If this solo particle can 
generate the required FCNC, then we are in Scenario B or Scenario A respectively. 



3.5 Scenario D: <y?p, ips ^ {0, tt} 
3.5.1 General Formulation 

In this scenario we assume no CP violating phases in the Bg — )■ yU+yU^ decay mode: 
V^P; V^5 ^ {0, vr} This is equivalent to all of the Wilson coefficients taking real values. 

■^In MFV this is the case. Namely C'g p/Cs,p ~ ms/mb- 
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Figure 4: Scenario D: ipp,ips £ {0,7r}. The correlation between the \P\ and \S\ param- 
eters for varying values of ^^r- ^^^^ shown is the current measurement of R. 



Clearly this constraint can also be applied to the other scenarios discussed in this section, 
but this scenario is distinct in that S and P are allowed to remain arbitrary real values. 
Yet, in the presence of a non- vanishing NP phase (f)f^, the CP-asymmetry S^^ could be 
non- vanishing. 

The resulting time dependent observables in this scenario are 



Ar 



cos 



-NP 



|P|2- 



|P|2 + |^|2 ' 

and the branching ratio observable is given by 

1 +ys cos 0f ^ 



— sm 



IPI 



151 



|P|2+ |5P 



(49) 



R 



+ \s\ 



I -Vs COS 



(50) 



Importantly, whereas the branching ratio observable R gives their squared sum, the 
is sensitive to the difference. With known these three observables are sufficient 
to determine the two NP unknowns \P\ and IS*]. As sin0^^ is already known to be 
small, Sfj_^ is also small in this scenario. Consequently A'^^ and R will be the relevant 
observables in this determination. With cos very close to unity one finds then 



IPI 



[l + Vs 



R 



Ar 



\S\' 



[l + Vs 



R 

'2" 



1 - ^^r 
l + VsAZ^^ 



(51) 



Finally a measurement of incompatible with the known value of would exclude 
this scenario and indicate new CP violating phases in the decay. 



In Figure 4 we illustrate how measurements of R and "^^^ be used to pinpoint 



the parameters IS"! and |P| (we have taken 



0). 
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3.5.2 Example of Models 



Minimal Flavour Violation (MFV) 

Models with Minimal Flavour Violation (MFV), but without flavour blind phases, as 
formulated as an effective field theory in Ref. |47J, belong naturally to this class. MFV 
protects against any additional flavour structure or CP violation beyond what is already 
present in the CKM matrix, while still allowing for additional, higher-dimensional, op- 



erators 47 . MFV therefore falls into Scenario D, with the added restriction that also 
is zero. Thus in models with MFV, as seen in (51), the time-dependent untagged 



observable together with the branching ratio observable R are sufficient to disen- 
tangle the scalar contribution S from P. A measurement of iS^^ ^ would falsify MFV. 
Typical examples in this class are MSSM with MFV and 2HDM with MFV. 

An exception are models with MFV and flavour-blind phases, like the 2HDM with 



such phases, also known as 2HDMj^jpY 48 . In this case model specific details are nec- 



essary in order for the time-dependent observables to distinguish between the operators 
and phases. 

3.6 Scenario E: P = 

In this scenario Cp\ C[q or 6Cio destructively interfere with CfJ^ to drive P to zero. 
Then non-zero values of the Bs — )■ ^~ observables will be driven purely by the opera- 
tors of. 

This scenario is similar to Scenario A, in that there is no rate asymmetry between 
the individual helicity decay modes. Thus the time-dependent observables are not inde- 
pendent: 

= - cos(2^s - <PT), = - sm(2^5 - <PT)- (52) 

The key difference, however, is that now only a scalar and not a gauge boson or pseu- 
doscalar is at work. Moreover, in the absence of new CP- violating phases A'^^ ~ 
which differs by sign from the Standard Model value and the analogous case in Sce- 



nario A as seen in (42). This is also a limiting case of Scenario D. The branching ratio 
observable is given by 

R=\S\^\ '-y^'''^^'''-^^"^ . (53) 

1 + Z/s 

We do not know any specific model that would naturally be placed in this scenario but 
we will investigate in Section |4] whether requiring tree-level exchanges of a pseudoscalar 
to cancel SM contribution is still consistent with the data. 



3.7 Summary 

In Table [T] we collect the properties of the selected models discussed above with respect 
to the basic phenomenological parameters listed in (36) and the class they belong to. We 
also indicate whether the phase can be non-zero in these models. In all cases \P\ is 
generally different from zero as it contains the SM contributions. In order to distinguish 
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Table 1: General structure of basic variables in different NP models. The last three cases 
apply also to the MSSM. 



between different models in each row of this table a more detailed analysis has to be 
performed taking all existing constraints into account. However, already identifying 
which of these four rows has been chosen by nature would be a tremendous step forward. 

4 Specific Models and Constraints from Bg Mixing 

4.1 Tree-Level Neutral Gauge Boson Exchange 
4.1.1 Basic Formulae 

As the first class of specific models we consider Z' models in which NP contributions 
to FCNC observables are dominated by tree-level Z' exchanges. A detailed analysis of 



these models has recently been presented in Ref. 14 . Also there the three observables 
in ^ have been considered but the emphasis has been put on the correlations of them 
with AF = 2 observables, in particular 5*^0. Here we will complement this study by 
computing the correlations among R, .4 at, and Snu, while taking the constraints from 



AF = 2 observables obtained in Ref. 14 into account 



We define the flavour-violating couplings of Z' to quarks as follows 

CfcMZ') = [At{Z'){s^,PLb) + At{Z'){s^,Pnb)] (54) 

where Af^p{Z') are generally complex. 

We also define the Z' couplings to muons 

C^Z') = [A'liZ')ih,Pd) + A'^iZ')ih,Pni)] Z'^ (55) 

and introduce 

Af{Z') = A^^(Z') - Af{Z'). (56) 

Then the non-vanishing Wilson coefficients contributing to Bg — )■ are given as 
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follows: 



where 



sin^ OwC 



-r]YYo{xt) - 



1 1 



Af{Z')A>f{Z') 



91m Ml VZV,, 



91m Ml v:y,. 



2 



a 



' y/2 2tt sin^ 6 



w 



(57) 
(58) 

(59) 



As only the coefficients Cio and C[q are non-vanishing this NP scenario is governed 



by the formulae (38 ) and (39 ). Indeed this scenario is an example of Scenario A in which, 
in addition to S' = 0, also the pseudoscalar contributions vanish. Yet, as P can differ 



from unity and have a nontrivial phase, a rich phenomenology is found 14 . 



4.1.2 Numerical Analysis 

It is not our goal to present a full-fledged numerical analysis of all correlations including 
present theoretical, parametric and experimental uncertainties as this would only wash 
out the effects we want to emphasize. Therefore we simply choose the three parameters 
entering our formulae, Fb,, t{Bs) and \Vts\ to be in the ballpark of their present central 
values: 

Fb, = 225.0 MeV, r(5,) = 1.503 ps, \Vts\ = 0.040. (60) 



Other relevant input can be found in the Tables of Ref. 31 . 

The main theoretical uncertainties in our analysis are due to the constraints on the 
couplings Af j^{Z) and Afj^{H) coming from the experimental values of AMg and 3^,^. 
Indeed the hadronic matrix elements of new operators are still subject to significant 
uncertainties. We will not recall the relevant formulae as they can be found in Ref. [14]. 
We will use the full machinery presented in that paper, setting the relevant parameters 
at their central values and requiring AMs and 5*^^ to be in the ranges 

16.9/ps < AM, < 18.7/ps, -0.20 < S^^ < 0.20. (61) 

Concerning the first range, it effectively takes the hadronic uncertainties into account. 
The second range corresponds to the 2a range for c/), in ([T]). 

As far as the direct lower bound on Mz' from collider experiments is concerned, the 
most stringent bounds are provided by CMS experiment |49] . The precise value depends 
on the model considered. While for the so-called sequential Z' the lower bound for Mz' 
is in the ballpark of 2.5 TeV, in other models values as low as 1 TeV are still possible. 
In order to cover large set of models, we will choose as our nominal value Mz' = 1 TeV. 
With the help of the formulae in Ref. fl^ it is possible to estimate approximately, how 
our results would change for 1 TeV < Mz' < 3 TeV. 



As in the analyses in Refs. 14,31 it will be instructive to consider the following four 



schemes for the gauge boson couplings and in the next subsection for scalar couplings: 
1. Left-handed Scheme (LHS) with complex A^'' ^ and A^^ = 0, 
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J.p I.U l.b 

R = BR(B, ^ ;i+;i-)/BRsii(B, ^ /i+fr) 



J.p I.U i.s 



Figure 5: iS^^j versus i? for LHS (left) and RHS (right), assuming Mz' = 1 TeV and 
A'^{Z') = 0.5. Gray region: exp la range for R. 



2. Right-handed Scheme (RHS) with complex A^^ ^ and A^^'' = 0, 

3. Left-Right symmetric Scheme (LRS) with complex A^f = A^ ^ 0, 

4. Left-Right asymmetric Scheme (ALRS) with complex A^f = — A^ 7^ 0. 

Note that the ordering in flavour indices in the couplings in these schemes is governed 
by the operator structure in mixing 14,31 and differs from the one in (54) and 

AfJH) = [A%,iH)Y, (62) 



(66). In this context one should recall that 



where H stands for either scalar or pseudoscalar. 



The ranges for Bg mixing given in (61 ) result in two allowed regions for the magnitudes 
and phases of the quark couplings A|^ depending on the scheme chosen above. These 
regions in parameter space are dubbed oases. The oases for each case have a two fold 
degeneracy in the complex phase of the coupling. Where it is relevant we will distinguish 
between these two different oases using the colours blue and red. 

In order to perform the present analysis we assign A^'^(Z') = 0.5, as was done 
in Ref. Il4j. In Section 4.2.3 and beyond, where we compare Z' exchange with various 



(pseu do) scalar exchanges, this coupling will be allowed to vary. The sign of this coupling 
is crucial for the identification of various enhancements and suppressions with respect 
to SM branching ratio and CP asymmetries and impacts the search for successful oases 
in the space of parameters that has been performed in Ref. 14 . If the sign of the Z' 



coupling to muons will be identified in the future to be different from the one assumed 
here, it will be straightforward, in combination with the discussion in Ref. 14 , to find 



MM 



out how our results will be modified. In Figure |5] we show the correlation between 
and R for LHS (left) and RHS (right). Corresponding correlations between and R 
and between A^^^ ^^"^ '^mm given in Figure [6]for LHS only. The two colours correspond 
to two oases in the values of the coupling Al^^IZ') that are consistent with AMg and 
5^0 constraints. 
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Figure 6: versus R (left) and v4^p versus S^^ (right) for LHS, assuming Mz' = 1 TeV 
and A^^(Z') = 0.5. Gray region in left panel: exp la range for R. 



We observe in analogy with findings of Ref. 14 that the correlations in the LHS 
and RHS schemes have the same shape except the oases and consequently the colours in 
Figure [5] have to be interchanged. We conclude therefore that on the basis of the three 
observables considered by us it is not possible to distinguish between LHS and RHS 
schemes because in the RHS scheme one can simply interchange the oases to obtain the 
same physical results as in LHS scheme. Consequently if one day we will have precise 
measurements of ^^p, S^+^i- and BR(i?s — )■ we will still not be able to distinguish 

for instance whether we deal with LHS scheme in the blue oasis or RHS scheme in the 
red oasis. 



As pointed out in Ref. 14 in order to make this distinction one has to consider 
simultaneously B — )■ K* ix~ and B — )■ Kfi~^fi~ decays and b — )■ suu transitions but this 
is beyond the scope of our paper. 

Yet, the plots in Figures [5] and [6] do give us valuable information. The allowed values 
for the three observables considered are, in this NP scenario, 

0.4 < AZ^ < 1.0, 0.2 < 1 < 0.9, 0.5 < i? < 1.8. (63) 

Moreover, the smallest values of A'^ and largest values of 1 are obtained for small- 

est values of R. The non-zero values of Sf^+^- originate in Z' models from requiring that 
AMg is suppressed with respect to its SM value in order to achieve a better agreement 
with data. As we will see below for models with scalar or pseudoscalar exchanges, this 
requirement can also be satisfied for a vanishing . 

If both LH and RH currents are present in NP contributions, but we impose a symme- 
try between LH and RH quark couplings, then NP contributions to Bs — )■ fi'^fi~ vanish. 
We thus find that 

^^'^r = cos(0f), 5,, = -sin(0f). (64) 
The branching ratio observable is given by 

"l + y,cos(0r) 



R 



(65) 
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which can also be obtained from (38) and (39) by setting P = 1. In view of the smallness 
of the results for the three observables are very close to the SM values. Still this 
example shows that even if no departures from SM expectation will be found in — )■ 
fi~^fi~ this does not necessarily mean that there is no NP around as with left-right 
symmetric couplings this physics cannot be seen in this decay except for small effects 
from the Bs mixing phase. This physics could then be seen in S — )■ K*fi~^fi~ and 
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B — 7- Kfi~^fi~ and b — )■ suu transitions as demonstrated in Ref. 

In the ALRS scheme NP contributions to Bg — )■ enter again with full power. 

Therefore the three observables in (|8| offer as in the LHS and RHS schemes good test of 



NP. In fact, as found in Ref. 14 , after the AB = 2 constraints are taken into account the 
pattern of NP contributions is similar to LHS scheme except that the effects are smaller 
because the relevant couplings have to be smaller in the presence of LR operators in 
AB = 2 in order to agree with the data on AMg. Therefore we will not show the plots 
corresponding to Figures |5] and [6j 

4.2 Tree-Level Neutral (Pseudo) Scalar Exchange 
4.2.1 Basic Formulae 

We will next consider tree-level pseudoscalar or scalar exchanges that one encounters in 
various models either at the fundamental level or in an effective theory. We will denote 
by H any spin particle, and will refer specifically to a scalar or pseudoscalar as or 
A^, respectively. It could in principle be the SM Higgs boson, but as the recent analysis 



in Ref. 31 shows, once the constraints from AF = 2 processes are taken into account, 
it is not possible to generate measurable effects in Bg — )■ through a tree-level 

SM Higgs exchange. The SM Higgs coupling to muons is simply too small. Therefore, 
what we have in mind here is a new heavy scalar or pseudoscalar boson encountered 
in 2HDM or supersymmetric models. Yet, in this subsection we will make the working 
assumption that either a neutral scalar or pseudoscalar tree-level exchange dominates 
NP contributions. A general analysis of FCNC processes within such scenarios has been 



recently presented in Ref. 31 . Also the observables in ([8j) have been analysed there, but 
with the emphasis put on their correlations with AF = 2 observables, in particular S^^. 
Here we will complement this study by computing the correlations among R, and 
iS^^, while taking the constraints from AF = 2 observables computed in Ref. 31 into 
account. 

We define the flavour violating couplings of H as follows 

CfcMH) = [Af{H){sPLb) + At{H){sPRb)] H (66) 
where Af^j^[H) are generally complex. Muon couplings A^'^^(if) are defined in a similar 



way. Note that through (62) in LHS and RHS schemes only A|j(if) and Af^{H) are 
non- vanishing, respectively. 
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Then the relevant non-vanishing Wilson coefficients are given as follows 
777,sin2^^C^ = -^T7^ VV ' 

where we have introduced 



Af{H) = A^^(iJ) - Af{H). 



(71) 



From the hermiticity of the relevant Hamiltonian one can show that A^'^(if) is real 
and A^p[H) purely imaginary. For convenience we define 

Af{H) = tAf{H), (72) 

so that A^[H) is real. 

Already at this stage it is instructive to see how different scenarios introduced in 
Section |3] are realized in this case. 

Scenario A: 

In this scenario 5 = 0. This can be realized simplest by setting A'g^{H) = 0, which, 



through (71), implies that A^{H) must be non- vanishing if the quark couplings differ 
from zero. The second possibility are left-right symmetric coupling to quarks but this 
would automatically imply also the vanishing of pseudoscalar couplings giving P = 1. 

Scenario B: 

We have just seen how this scenario can be obtained as a limiting case of scenario A. In 
order to have non- vanishing S in this case this is realized by setting A^p{H) = 0, which 



through (71) implies that A^'^(if) must be non- vanishing if the quark couplings differ 
from zero. 

Scenario C: 

To achieve this scenario the scalar coefficients should be equal up to a sign to the 
pseudoscalar ones. This requires the exchanged spin-0 particle to be a mixed scalar- 
pseudoscalar state, which is beyond the scope of the present analysis. We will instead 
realise Scenario C in Section 4.3| by considering the presence of both a scalar and a 



pseudoscalar with equal masses and equal couplings to quarks. 
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Figure 7: The dependence of the observables R (left) and A^^y (right) on the pseudoscalar 
lepton couphng Ap^{H) satisfying the Bs mixing constraints in the LHS case. For a 
pseudoscalar with mass M^o = 1 TeV. 

Scenario D: 

Because a single scalar or pseudoscalar allows only 5 or P to deviate from its SM value, 
respectively, the intended usage case of this scenario, namely arbitrary but real valued 
S and P, cannot be realised. 



Scenario E: 



(') 

In this concrete model in which there are no NP contributions to C^q the vanishing of 
P implies: 



m 



Bs 



2m, 



rrib 



C'p — Cp 

•-^10 



1. 



(73) 



We will investigate whether this condition is consistent with existing constraints when 



the relevant Wilson coefficients are given as in (69) and (70) 



4.2.2 Numerical Analysis 



Analogous to the case of tree-level Z' exchanges we will use the results of the AF = 2 
analysis in Ref. |31| to constrain the quark-scalar couplings in the schemes LHS, RHS, 

The next step is to set values 
For a single scalar particle the 



LRS and ALRS by imposing the conditions in (61). 
for the scalar and pseudoscalar muon couplings. 



parameter \S\ driving NP {Scenario B) is directly proportional to the muon coupling 
\A^g'{H^)\. However, for a single pseudoscalar particle ^4°, the muon coupling Ap^(A°) 
is not directly proportional to P, and the resulting NP observables thereby have a more 
involved dependence on it. In Figure M we show the dependence of the observables R 



(left panel) and v4^p (right panel) with respect to muon coupling A^p{A^) defined in (72) 
satisfying the B^ mixing constraints for the LHS case. We observe that the parameter 
space of the NP physics observables is very dependent on whether we pick a large or 
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Figure 8: iS^^ versus R for for LHS scheme with a scalar (left) and pseudoscalar (right) 
for Mho = Mao = 1 TeV. Gray region: exp la range for R. 



small coupling, and that a fixed coupling cannot do it justice. We further observe that 
the oases become indistinguishable if the sign of the coupling is not fixed. 

In order to compare the oases behaviour of the scalar and pseudoscalar we begin by 
fixing the muon couplings to 

Af (i7°) = 0.040, Af (A°) = i 0.020, (74) 

and AfiH'^) = Af (A^) = 0. 



As demonstrated in Ref. 31 these values are consistent with the allowed range for 



B{Bs — )■ fi~^fi ) when the constraints on the quark couplings from — are taken 



into account and M = 1 TeV. All other input parameters are as in Ref. 31 . The reason 
for choosing the scalar couplings to be larger than the pseudoscalar ones is that they are 
more weakly constrained than the latter because the scalar contributions do not interfere 
with SM contributions. 

In Figure sjwe show the correlations of versus R satisfying Bg mixing constraints 
for a single tree-level scalar (left) and pseudoscalar (right) exchange in the LHS scheme. 
For the scalar case the blue and red oases overlap. The red oases in the pseudoscalar 
case corresponds to i? < 1 and is therefore clearly distinguishable from the scalar case. 



where R > 1 for both oases. In Section 4.2.3 we will compare these correlation with Z' 
exchange. 

As we stated earlier, fixing the pseudoscalar muon couplings to one value does not 
reveal the full structure of the NP parameter space. We therefore now consider the muon 
couplings varied over the following range: 

\Af{H% |Af (A°) I G [0.02,0.04]. (75) 

From here on we will ignore the sign of the lepton couplings, but again note that this 
degeneracy can be resolved in the pseudoscalar case if the blue or red oasis from the Bs 
mixing constraints can be singled out. We thus also stop distinguishing between the two 
oases. 
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Figure 9: Overlay of the correlations for R versus (left) and S^^ (right) for 

the various specific models considered. The lepton couplings are varied in the ranges 
\A''/p{H)\ e [0.02,0.04] and A''/{Z') E [0.3,0.7]. All particles are taken to have a mass 
of ITeV. 



In Figure [q , i? is plotted against (left panel) and (right panel) with the 
regions allowed by the Bg mixing constraints overlayed for the various specific tree- 
level models discussed in this section. The scalar and pseudoscalar muon couplings 
have been varied as just discussed, and also the Z' muon couplings have been varied: 
A'^{Z') G [0.3, 0.7]. These three models are shown for the LHS scheme. Focusing for the 
moment on models with single scalar (H^) or pseudoscalar (A^) exchanges, the following 
observations can be made: 

• The branching ratio observable R can in the scalar case only be enhanced as there 
is no interference with the SM contribution. On the other hand, in the pseudoscalar 
case it can be suppressed or enhanced depending on which oasis of parameters is 
chosen. 

• The values of ^^p are positive for both and and for R within one a exper- 
imental value close to unity. 

• can reach ±0.50 in both cases. 

We do not show the corresponding results in the RHS scheme as, similarly to the 
gauge boson case, the correlations in question have identical structure with the following 
difference between scalar and pseudoscalar cases originating in the absence and presence 
of correlation with SM contributions, respectively: 

• In the scalar case the two correlations in question are invariant under the change 
of LHS to RHS. 

• In the pseudoscalar case the structure of two correlations remain but going from 
LHS to RHS the colours have to be interchanged as was the case for gauge bosons. 
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In the LRS expected, NP effects are very small as scalar and pseudoscalar 



in (74): 



contributions are absent and (|64|) applies. We then find for the muon couplings fixed as 

(76) 



0.984 < ^^'^p < 1.00, 



< 0.18. 



Finally we investigated whether the relation (73), representing Scenario E is still 



consistent with all available constraints. This is not the case if we take the pseudoscalar 



lepton coupling chosen in (74) and a mass for the pseudoscalar of 1 TeV. For the LHS and 
RHS schemes a lepton coupling of Ap^(if) k, ±i 0.1 is needed to satisfy the relation. If 
a pseudoscalar does manage to make P vanish, then a scalar particle is needed to satisfy 
the lower bound on R. Such a model, with both a pseudoscalar and scalar particle 



present, is discussed in Section 4.3 



4.2.3 Comparison with Z' Scenario 

While the discussion presented above shows that the contributions of scalars and pseu- 
doscalars can be distinguished through the observables considered, more spectacular 
differences occur when one includes the Z' scenario in this discussion. Indeed the cor- 
relation between iS^^ and R in the left panel of Figure s] has a very different structure 
from the case of pseudoscalar or scalar exchanges shown in Figure |8j 

In the right panel of Figure |9] an overlay of these regions is shown for LHS schemes. 



with the lepton couplings varied as given in (75). Similarly, in the left panel of Figure |9] 
we show the correlation between ^^'^ and R, where strong contrasts between the allowed 
regions also emerge. The difference between the Z' and pseudoscalar exchange is striking 
because, unlike for a scalar, both particles generate Scenario A. 

The difference between the A^-scenario and Z'-scenario in question can be traced 
back to the difference between the phase of the NP correction to P, which was defined 



in (40). As the phase ^23 in the quark coupling A^'^ from the analysis of S^-mixing in 



Ref. 31 is the same in both scenarios the difference enters through the muon couplings. 



which are imaginary in the case of the A^-scenario but real in the case of Z' . This is 
why the structure of correlations in both scenarios is so different. Taking also the sign 
difference between Z' and pseudoscalar contributions to the h — )■ sfi~^fi~ amplitude into 
account, we find that 



P{Z') = 1 + rz'c' 



P(A°) = 1 + rAoe'- 



with 



rz' ~ rAo, 



523 - 



5ao = 82 



'aO 
TT 

2' 



(77) 



(78) 



It turns out that the Bg mixing constraints force the phase to be in the ballpark of 
90° and 270° for the blue and red oasis, respectively 14,31 . This implies in the case 



of the Z' scenario, as seen in Figure [sj positive and negative value of S^^ for the blue 
and red oasis, respectively. Simultaneously R, where NP is governed by cos6z', can be 



enhanced or suppressed in each oasis. On the other hand, (78) implies that the phase 



Sao is in the ballpark of 0° and 180° for the blue and red oasis, respectively. Therefore 
the asymmetry S^^ can vanish in both oases, while this was not possible in the Z' case. 
As NP in R is governed by cos(5^o, this enhances and suppresses R for blue and red 
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oasis, respectively as clearly seen in Figure [8j In particular, R differs from its SM value, 
while this is not the case in the Z' scenario. Finally, let us note that with larger values 
of muon couplings NP effects in R, and iS^^ can be larger than shown in Figs, [s 
and|6| (see, for example. Figure 10). 



What is particularly interesting is that these differences are directly related to the 
difference in the fundamental properties of the particles involved: their spin and CP- 
parity. As far as the last property is concerned, also differences between the implications 
of the pseudoscalar and scalar exchanges have been identified as discussed in detail above. 
They are related to the fact that the scalar contribution, being CP even, cannot interfere 
with the SM contribution. 



4.3 Tree-Level Neutral Scalar+Pseudoscalar Exchange 
4.3.1 Basic Formulae 

In this model we assume the presence of a scalar and pseudoscalar with equal (or 
nearly degenerate) mass Mh- This is, for example, effectively realised in 2IIDMs in a 
decoupling regime, where and are much heavier than the SM Higgs and almost 
degenerate in mass [43]. We will show that under specific assumptions this setup can 
reproduce Scenarios C, D or E. 

The couplings of the scalar and pseudoscalar to quarks are given in general by the 
following flavour-violating Lagrangian: 

£fcnc(^°,A°) = [Al\H'){sPLb) + At{H'){sPRb)] 

+ [At{A'){sPLb) + At{A'){sPnb)] A'. (79) 

We will assume that the scalar and pseudoscalar couple with equal strength to quarks: 

C 3 Dl ADn{H^ + iA'^)+ h.c, (80) 

where D = [d, s, b) and A is a matrix in flavour space. Then 

A'^(H^) = A'\ Af{H^) = 



.bs 



A^^(A°) = iA'\ Af{A^) = -i A"' . (81) 



, bs 



where in general A*^, A''* G C. 



Scenario C: 



To reproduce this scenario we set the pseudoscalar and scalar masses to be exactly equal: 
Mffo = Mao = Mh- Further relating the lepton couplings by a single real parameter 

A^^(//°) = A'^'', A^''(A°) =iA>'f' (82) 
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and inserting the lepton and quark couplings into formulae (67)-(70) we find: 

1 ^sb^fifj, 



Cq — —Cf 



L/C — L/ p 



A 



hs 



gl MMlmbSui^ew Vtl Vtb 



(83) 



(84) 



This simple model satisfies the relations in (48) and thereby belongs to Scenario C. 
These relations are in fact valid for all the quark coupling schemes: LHS, RHS, LRS and 
ALRS. Yet the physics implications depend on the scheme considered: 

• In LHS and RHS schemes NP contributions to B^-B^ mixing from scalar and 
pseudoscalar with the same mass cancel each other so that there is no constraint 
from Bg-Bg mixing. Thus NP effects in Bg — t- fi^ii~ can only be constrained by 
the decay itself or other b — )■ si'^£~ transitions. 



In LRS and ALRS schemes non-vanishing contributions from LR operators to B 
Bg mixing are present. Moreover we find 



Cs = C's, Cp = -C'p (LRS), (85) 

Cs = -C's, Cp = C'p (ALRS). (86) 

Therefore in the LRS case only pseudoscalar contributes to Bg — > fj^'^fi" {Scenario 
A), while in the ALRS case only scalar contributes {Scenario B). 



We conclude therefore that in order to have an example of Scenario C that differs 
from Scenario A and B and moreover in which NP contributions to B^-B^ mixing are 
present, we need both L and R couplings which are not equal to each other or do not 
differ only by a sign. 

An option to reproduce Scenario C with non-trivial constraints from mixing is given 
by Minimal Flavour Violation (MFV). In the MFV formalism A is constructed out of the 
spurion matrices Yjy and Yd [47] . In principle the following constructions can contribute 
to the 5 — )■ s FCNCs at leading order in the off-diagonal structure: 



u^u^D, YdYd^uYijYd, YjjY^YuY^Y]:). (87) 
However, the last two will in general receive dynamical (loop) suppressions. Thus, for 



YttYX Yr)Yn Yi 



simplicity, we assume the first construction to be dominant. In the notation of Ref. 48 



where MFV is discussed in the context of a general 2HDM with flavour blind phases 
(2HDMj,^jPY), this is equivalent to assuming |ao| 3> |ai|, |a2|. As a result we find 



tb, 



.bs 



vs v*v, 



tb 



^^_^sb 

rrib e 



(88) 



Thus under the above assumptions all of the quark couplings in (79) can be expressed 
in terms of a single NP parameter e |^ 



""^It should be emphasized that in general this is not the case for 2HDM j^py 
comments below. 
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See additional 
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The parameter e is real in pure MFV but may be complex in 2HDMj^jpY 48 . Inserting 



relation (88) into (84) we find 



C's='^-Cs C'p = -'^-Cp, (89) 

TUh e mb e 

and observe a nis/mh suppression of the primed operators. In pure MFV, where e is real, 
the parameters C^'p are also all real. 

Scenario D: 

In Scenario D the parameters P and S are arbitrary but do not carry new CP violating 
phases. The pure MFV model with a scalar and pseudoscalar that we just discussed 
is therefore a natural candidate. However, because this model was defined to satisfy 
Scenario C, as it stands we have P ± S = 1. If we continue to insist that the scalar and 



pseudoscalar should couple with equal strengths and phases to quarks as in (80), then 
there are two choices for making P and 5* arbitrary. 

One choice is to allow different couplings to leptons for the scalar and pseudoscalar i.e. 
\A'^P'{H^)\ 7^ |A'^'^(A°)|. In this case the constraints from Bg mixing (discussed below) 
do not change, and only the current bounds on R must be satisfied. 

Alternatively, a non-trivial difference between the scalar mass Mho and the pseu- 
doscalar mass Mao can be introduced. In this case the lepton couplings can remain 



equal as defined in (82). The catch, however, is that now the LL and (to a much lesser 
extent in MFV) RR contributions to Bg mixing no longer vanish. Thus the allowed mass 
difference, and thereby the arbitrariness of P and S is constrained by mixing. 

Scenario E: 

This scenario requires that P = and therefore that S alone generates a value of R large 
enough to meet the current experimental bounds. As we are dealing with two spin-0 



particles, the pseudoscalar in the present model must satisfy the relation given in (73). 



By definition Scenario E allows S to have a new CP violating phase (ps- However, 



the relation in (73) requires that arg(Cp — Cp) = 0. Therefore, because we required that 



the scalar and pseudoscalar should couple to quarks with equal strengths and phases (see 



(80)), it follows that (ps = 0. The Scenario E realisable in this model is therefore just 
a specific case of Scenario D, where P and S are real and arbitrary, with the addition 
that P is tuned to vanish. 

In the next section we will address whether, given that 5" and P are made to vary due 
to a scalar-pseudoscalar mass difference and P is tuned to zero, the range of S allowed 
by mixing can satisfy the experimental bounds on R. 

4.3.2 Numerical Analysis 

Our numerical analysis for this model will focus on the above mentioned assumptions 
that produce Scenario C. Specifically, we begin by assuming an exactly degenerate scalar 
mass Mfj, equal scalar and pseudoscalar lepton couplings and MFV. At the end of this 
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section we also briefly address the consequences of a scalar-pseudoscalar mass difference, 
which could produce Scenarios D and E. 



By imposing MFV on the flavour matrix A introduced in (80), it follows that the 
analogues of A** in the Ba and K systems are related to it by 



A 



db 



V* ~ 

td ^ sb 



A 



ds 



td 



A 



sb 



(90) 



Therefore the value taken by A*^ should in principle not only satisfy the experimental 
Bg mixing constraints, but also those of the Bd and K systems. In practice, however, NP 
contributions in this model to Bd mixing are suppressed by a factor of md/rris relative 
to Bg mixing and thereby very small. As a result, this model with MFV cannot relieve 
the current tensions in Bd mixing between theory and experiment 51,52 . Contributions 
to neutral Kaon mixing are totally negligible. We therefore proceed to only consider 
constraints from Bg mixing. 

The only contribution that survives in Bg mixing is the LR one and this introduces 
the following shift in the SM box function |31 



where 



[AS{Bg)],^ = 2r 



SiBg) = Soixt) + [ASiBg)]^j,, 

^^sb^^bs 



LR 



Mh{VtgV*y 



2r 



LR 



rUb 



Ml ■ 



(91) 



(92) 



with r^^ = -3 X 10^ TeV^ [31j. 

The following observations should be made: 



In spite of possible new flavour blind phases in the MFV scenario, these phases do 
not show up in i^^-niixing, so that the CP asymmetry S^^ remains at its SM value, 
still consistent with experiment. In Ref. 48 the Ai? = 2 operator for a 2HDM 
with MFV is also found to leave flavour-blind phases unconstrained in the limit 
|ao| ^ |ai|, |a2|. In general this is not the case for 2HDMj^jpY and, as analysed in 
Refs. (48l|50l 



5*^0 can receive NP contributions. Note that these phases also appear 
in the flavour-conserving Yukawa couplings, which contribute to the electric dipole 
moments of various atoms and hadrons by the exchange of Higgs fields. But as 

in Ref. 



shown for 2HDM 



MFV 



50 the present upper bounds on EDMs do not yet 



have any impact on the observables considered here 

On the contrary, AM^ receives a small (suppressed by m^/m^,) negative contribution 
which is good as the SM value is roughly 10% above its experimental value |lj . This 
suppression due to LR operators within a MFV framework was first pointed out 
for the MSSM with MFV in Ref. 
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The fact that the flavour-blind phases are unconstrained through i?s-niixing allows 
us to obtain significant effects from them in Bg — )■ observables as we will see 

soon. 



^We thank Minoru Nagai for enlighting comments on these issues. 
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Figure 10: Overlay of the correlations for R versus (left) and S^^ (right) for 

the various specific models considered. The lepton couplings are varied in the ranges 
\A''/p{H)\ e [0.00,0.06] and A''/{Z') E [0.0, 1.0]. All particles are taken to have a mass 
of ITeV. 



For both MFV and MFV we find the range: 

|A'^| G [0.00196,0.00530], 



(93) 



for Mh = 1 TeV, which, as seen in (88), is consistent with the tacit assumption that e 
should be small. 

To proceed with numerics for the — )■ fi^fi~ observables we must set the coupling 
of and A° to muons. In the context of Scenario C this means setting A'^'^ as defined 



in (82). In order to compare with the single tree-level scalar and pseudoscalar models 



discussed in the previous section we begin by varying the coupling between [0.02,0.04] 



as also done in (75) 



The left panel of Figure [9j shows A'^^' plotted versus R for MFV with = 1 TeV. 
The allowed region from Bg mixing constraints shown in this plot should be compared 
with the theoretical situation sketched for Scenario C in the left panel of Figure [3] 
By inspection of the theoretical plot one observes that the pure MFV model (with no 
flavour blind phases) corresponds to the outer border of the MFV region shown. It is 
interesting to observe that in both models a negative is possible within the Bs 
constraints mixing, in contrast to the tree-level models considered above with a single 
(pseu do) scalar or gauge boson. Because the flavour-blind phase in MFV is completely 
unconstrained, almost the entire experimentally allowed region is left unconstrained by 
Bg mixing in this model. 

In the right panel of Figure [9 we similarly show versus R in the MFV model 
for Mfj = 1 TeV. In the pure MFV model S^^ = and therefore these plots are not 
interesting. 

In a 2HDM with large tan /3, which can generate a decoupled heavy scalar and pseu- 
doscalar as discussed here, the muon coupling is given by 



'V2 



m 



^tan/3 



j = -0.05 


tan P 
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(94) 
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Mh [GeV] Mh [GoV] 



Figure 11: The allowed region of versus the heavy scalar mass Mh in MFV (left 
panel) and MFV (right panel). The allowed region satisfies the Bg mixing constraints 
and falls with the 2 a C.L region oiR: Re [0.30, 1.80]. 



Model 


MHo[GeV] 


S 


R 


MFV 
LHS, RHS 


900 - 1095 
1030 - 1310 


1.20 - 0.81 
0.94 - 0.56 


1.22 - 0.58 
0.26 - 0.75 



Table 2: Ranges allowed by Bs mixing for the Scalar+Pseudoscalar model with non- 
degenerate masses. The model has been tuned so that P = 0, and thereby corresponds 
to Scenario E. The pseudoscalar mass is Mao = 1 TeV. 



which demonstrates that the (pseudo)scalar muon couplings can be larger than what we 



have assumed so far. In Figure 10 we repeat the plots we have shown in Figure p\ but 



now with the muon couplings varied over much larger ranges: 

|Ag|'p(i7°, A°) I e [0.00,0.06], |A^^(Z')| G [0.0,1.0]. (95) 

This range of couplings gives a better impression of the full allowed parameter space, at 
the cost of hiding some of the characteristic differences between the considered models. 
We do not again show the large allowed region of H° + A° model with MFV, but we do 
show it with pure MFV in the A'^^' versus R case (left panel). 



In Figure 11 we show the allowed range of A'^^ with respect to the heavy scalar mass 



Mh in the MFV (left panel) and MFV (right panel). In these plots we have fixed the 
muon couplings to A^'^ = —0.05. The allowed range takes the Bg mixing constraints 



into account and falls within the 2 a C.L of R as defined in (29). We observe that 
for Mh < 0.75 TeV negative values of A'^y' predicted in this scenario, while for 
Mh > 2.5 TeV its SM value is approached. 

We will now consider a mass difference between the scalar and pseudoscalar: 



Mho ^ Mao. (96) 

We consequently shift our focus from Scenario C to Scenario E. Note, however, that 
a small mass difference is still consistent with a 2HDM in a decoupling regime, and 
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will approximately reproduce Scenario C. In the presence of a non-zero mass difference 
the following contribution to the SM box function in MFV or a LHS model no longer 
vanishes: 



[AS{B,)] 



[A 



sb~\ 



LL 



HO] 



XL 



(M 



AO. 



(97) 



and an analogous expression for a RHS model after the replacement [A**^]* — )■ A*'' and 
L with ~R with r^^ = r^^ = 50 TeV^ In MFV we therefore find: 

We observe that: 



rrib 



i[^SiB,)Ur = {erylyt 



LL/ 



(98) 



The LL contribution, while suppressed relatively to the LR contribution through 
smaller hadronic matrix element (|r-^-^| ^ 6|r-^-^|) and the splitting between scalar 
and pseudoscalar masses, does not suffer from m^/mb suppression. Thus whether 
the SLL contribution or LR dominates depends sensitively on the size of the mass 
splitting in question. The SRR contribution is totally negligible. 

The SLL contribution now contains in principle a new flavour-blind phase in e 
allowing for new CP-violating effects in Bg mixing. 



The formulae (92) and (98) are also valid in a non-MFV framework in which new 



flavour and CP-violating phases are present in e. 
For our numerics we flx the pseudoscalar mass to M^o = 1 TeV and will continue 



to assume a universal lepton coupling as given in ( 94 ) . The quark coupling required to 
set P = has strength |A"^| = 0.0043 and a phase arg(A"^) = 0^^/2 = -1°. This 
then poses the question of what the allowed range for the scalar mass M^o is that is 
compatible with the Bg mixing constraints, and whether the resulting S satisfles the 
experimental bounds on R. 

Because the scalar and pseudoscalar couple in the same way to quarks (see (80)), 
P = implies that there are no new CP violating phases in the decay. The NP mixing 
phase (j)^^ can still contribute (unless we assume MFV), and we flnd for the time- 
dependent observables: 



Ar 



— cos 



iNP\ 



MM 



sm 



-NP^ 



(99) 



This change in sign for these observables with respect to the SM is a smoking gun signal 
of scalars dominating in the Bg — )■ fi^fj,^ decay. In Table [2] we summarise our results for 
the allowed ranges of the scalar mass M^o, parameter 5* and observable R for both MFV 
and a LHS /RHS quark model. Both models are seen to satisfy the current experimental 



range for R given in ( 29 ) . 
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5 Summary 



We have performed the first detailed phenomenological analysis of the time-dependent 
rate for the Bg — decay following the formalism developed in Ref. [s]. Our anal- 
ysis demonstrates that decay-time studies of Bs — )■ Ai"'"yU~, which offer the observables 
and Sfj^n in addition to the branching ratio, allow for various NP scenarios to be 
disentangled. Specifically, the presence of new scalar, pseudoscalar or gauge boson par- 
ticles can potentially be identified, which is not possible on the basis of the branching 
ratio alone. 

We have proposed a classification of various NP scenarios in terms of two complex 
variables S and P that fully describe the three observables involved, and can be ex- 
pressed in terms of the fundamental parameters of a given model. The experimental 
determination of S and P, accompanied by plausible model specific assumptions, will 
allow us to probe NP in this theoretically clean decay. We have illustrated this by placing 
several popular extensions of the SM into phenomenological scenarios introduced by us 
see Table ig. 

We have further presented numerical analyses for the observables in question in mod- 
els for tree-level contributions to Bg — t- mediated by heavy gauge bosons, scalars 
and pseudoscalars. The plots in Figures [9] and 10 illustrate our general findings. Our 
main messages from these analyses are as follows. 

• The phenomenology of a tree-level Z' exchange with respect to the studied ob- 
servables is very different in structure to that of spin-0 particles, in particular 
pseudoscalars. This is shown in Figure [9] (see also Figure IS] versus Figure Isl). 



• In turn, the phenomenology of a scalar is more restricted than that of a pseu- 
doscalar. For instance, suppression of R with respect to its SM value would exclude 
a NP scenario with only a single scalar, whereas such a suppression is possible for 
a single pseudoscalar. 

• For models with a single new particle with a mass of 1 TeV - specifically a gauge 
boson, scalar or pseudoscalar - negative values of require large couplings to 
muons and a significant deviation of the Bg mixing phase (pg from its SM value. 

• On the contrary, a negative value of can naturally be explained in models with 
both a scalar and pseudoscalar and a common mass Mh < 1.5 TeV. An example 
is decoupled 2HDMs. In this setup a deviation of (j)g from its SM value is not 
required. Furthermore, in these models R has a strict lower bound. 

• For a Z' scenario the required suppression of AMg, due to current tensions with 
experiment, implies that 7^ 0. 

• For a single pseudoscalar or scalar scenario, the required suppression of AMg im- 
plies a departure from the SM value of BR{Bg — )■ /i^/i"). 

The numerous plots and examples presented by us provide a roadmap for future 
experimental results of this outstanding rare B decay. 
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